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EFFECT OF PLUG AND SHROUD GEOMETRY VARIABLES ON 
PLUG-NOZZLE PERFORMANCE AT TRANSONIC SPEEDS* 

By Bobby L. Berrier 
Langley Research Center 

SUMMARY 

An investigation of the effect of plug and shroud geometry variables on the perfor- 
mance characteristics of a basic 15° cone plug nozzle has been made in a transonic wind 
tunnel. Geometry variables were plug truncation (0, 30, 72.5, and 100 percent of full 
length from annular throat), shroud boattail angle (8°, 20°, and 30°), and truncated plug 
base shape (flat and semitoroidal). All nozzle configurations were designed to operate 
at a jet total-pressure ratio of 16.5. The ratio of annular throat area to maximum cross- 
sectional area was 0.25. The wind-tunnel Mach number was varied from 0 to 1.30. The 
jet total -pres sure ratio was varied from 1.0 to approximately 8.0 depending on the exter- 
nal airstream Mach number. Model angle of attack was 0° for all tests. A comparison 
of experimental plug pressures and theoretical plug pressures obtained by the method of 
characteristics is made. 

The performance of the 15° cone plug nozzle was shown to be higher than the per- 
formance of a comparable isentropic plug nozzle. The results also show that plug trun- 
cation decreases cone-plug -nozzle performance although the 15° cone plug was truncated 
30 percent with a loss of generally less than 1 percent. The 72.5- and 100-percent- 
truncated plug-nozzle configurations show severe losses in nozzle performance. 
Increasing boattail angle generally decreased cone-plug-nozzle performance except for 
the 100-percent -truncated plug-nozzle configurations. 

INTRODUCTION 

As a result of relatively low lift-drag ratio and high net specific fuel consumption, 
airplanes which operate in the supersonic speed range require exhaust nozzles which 
operate at near optimum efficiency in the supersonic cruise condition. The payload of 
an airplane cruising at Mach 2.2, for example, is decreased by about 8 percent for a 
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loss of 1 percent in nozzle efficiency. (See ref. 1.) In addition, high nozzle efficiency is 
desirable for acceleration through the transonic speed range where airplane drag is high 
and during the cruise and loiter conditions at subsonic speeds. 

The wide range of flight speed and jet total -pressure ratio (i.e., the ratio of jet- 
exhaust total pressure to external stream static pressure) required for supersonic air- 
planes thus creates a serious problem in nozzle development for such airplanes. The 
simple convergent-divergent nozzles with fixed geometry give high nozzle efficiency 
when operating at the design pressure ratio. However, a decrease or increase of jet 
total -pressure ratio from the design jet total-pressure ratio causes a rather large loss 
in nozzle efficiency due to overexpansion or underexpansion losses. Further, at forward 
speeds, the airplane nozzle system incurs an additional loss due to pressure drag on the 
boattail surfaces. Overexpansion and underexpansion losses may be eliminated by pro- 
viding variable nozzle geometry through mechanical means. This solution leads, how- 
ever, to a complex design with fabrication difficulties, an increase in weight, and an 
increase in boattail pressure drag at low jet total -pressure ratios. 

Hence, lightweight exhaust nozzles that can provide high nozzle performance levels 
over a wide range of operating pressure ratios are needed. One promising nozzle type 
is the plug nozzle for which the outer boundary of the exhaust stream continually adjusts 
to ambient or free-stream pressure. 

The objective of this investigation was to determine the effects of plug and shroud 
geometry variables on the performance of a 15° cone plug nozzle operating at off -design 
jet total-pressure ratios in the subsonic and transonic speed regimes. Geometry vari- 
ables were plug truncation (0, 30, 72.5, and 100 percent of full length from annular 
throat), shroud boattail angle (8°, 20°, and 30°), and truncated plug base shape (flat and 
semitoroidal) . In this investigation, only plug nozzles with all exhaust flow expansion 
occurring external to the nozzle shroud were considered. 

Theoretical pressure distributions on the plug surface, as calculated by the method 
of characteristics, are compared with experimental pressure distributions at an external 
airstream Mach number of 0 and for various values of jet total-pressure ratio. 

The experimental investigation was conducted in the Langley 16 -foot transonic tun- 
nel at static conditions and at Mach numbers of 0.50 to 1.30 at an angle of attack of 0°. 

The jet total-pressure ratio was varied from 1.0 (jet off) to about 8.0, depending on the 
Mach number. The nozzles were designed to operate at a jet total-pressure ratio of 16.5. 
The jet exhaust was simulated by the decomposition products of 90-percent-concentration 
hydrogen peroxide which had a specific -heat ratio of 1.266 at a temperature of 1013° K. 
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SYMBOLS 


'D,/3 


cross-sectional area, meters^ 

b 

'C' 

boattail pressure drag coefficient > -Cp 


where b = -1.113d 


max 


H 

- 0 max 

for afterbody A 


b = -0.550d max for afterbody B 


b = -0.343d max for afterbody C 


C F,i 


thrust coefficient, 


q A. 


oomax 


F i F i 

ideal thrust coefficient, — - ; for static conditions, — 

l^max Poo max 


'(F-D) 


F D F-D 

thrust-minus-drag coefficient, - — ; for static conditions, 


q A, 


oo max 


Poo A max 


%P 


1.0 


C r -i plug thrust coefficient, 7 — 2-i. 

F >P lu g ' L A m!lv 

x A 1 = o 


max 


. x 11 xi.:, Boattail skin -friction drag 

skin -friction drag coefficient on nozzle boattail, c - 

1co A max 


Cf C yi skin -friction drag coefficient on cylindrical part of afterbody, 

Cyli nder skin -friction drag 
tLoo^-max 


P; " Poo 

pressure coefficient, -L- ; for static conditions, 


Pf -P. 


D boattail pressure drag and skin-friction drag, ^C D ^ + C f ^q eo A max , newtons 
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diameter, meters 


jet thrust, newtons 

general length or reference length of plug from shroud lip, positive down- 
stream, meters 

length of plug from throat to plug tip, 19.51 centimeters 
amount of plug length cut off, meters 
free -stream Mach number 

jet -exit exhaust Mach number for fully expanded flow 

measured mass flow rate, kilograms/second 

static pressure, newtons/meter2 

total pressure, newtons/meter^ 

dynamic pressure, newtons/meter^ 

gas constant, joules/ (kilogram) (degree Kelvin) 

radius, meters 

radius on plug base to pressure orifice location, meters 
radius defining semitoroidal-base plug geometry, meters 
temperature, degrees Kelvin 

axial coordinate from shroud lip, positive downstream, meters 

axial coordinate from shroud attachment point ^x/d max = -0.7 j, positive 
downstream, meters 


axial coordinate from plug attachment point, meters 



radial coordinate, meters 


external boattail lip angle, degrees 


ratio of specific heats 


internal boattail angle, degrees 
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Subscripts: 


flow deflection angle, degrees 


Mach angle, sin p = — , degrees 


b 

bal 


max 


0 


radial angle of nozzle, degrees 

afterbody 
base region 
balance reading 
nozzle exit 
internal or ideal 
jet flow 
local 

maximum 

condition at surface 
stagnation 
boattail 
free stream 


Superscript: 


* nozzle throat 

A bar over a symbol indicates an average condition. 

EXPERIMENTAL ARRANGEMENT 
Tunnel 

The experimental investigation was conducted in the Langley 16-foot transonic 
tunnel. The facility is a single -return, atmospheric tunnel with air exchange for cooling 
and has a continuously variable speed range from a Mach number of 0.20 to 1.30. The 
average Reynolds number per meter varies from 9.71 x 10® at a Mach number of 0.50 
to 12.6 x 10® at a Mach number of 1.30. Additional information about the facility is given 
in references 2 and 3. 


Model and Support System 

The exhaust stream of a turbojet engine was simulated by a hydrogen peroxide mono- 
propellant rocket system similar to the jet simulator described in reference 4. A sketch 
illustrating the major components of the jet simulator is given in figure 1. 

The strut -supported, turbojet -engine simulator model used in this investigation is 
illustrated in the sketch in figure 2, and a photograph of the model installed in the test 
section is shown as figure 3. The model consisted of a conical forebody, a cylindrical 
section of 15.24-cm diameter, and an afterbody -plug combination consisting of a cylindri- 
cal section, a boattail, and a cone plug. The afterbody -plug combinations were detach- 
able at the 104.39-cm station to facilitate configuration changes. All external surfaces 
downstream of this station were attached to the thrust balance. All the plug-nozzle con- 
figurations investigated were designed, by using the specific -heat ratio for the exhaust 
gases (y = 1.266), for a jet total -pressure ratio of 16.5 (M e = 2.46). A streamline is 
assumed to be parallel to the axis of symmetry and to be issuing from the boattail lip at 
design conditions. The ratio of throat critical area to maximum cross-sectional area 
was designed to represent a nonaugmentation condition (no afterburning). 

Eleven configurations were investigated by interchanging different shrouds and 
plugs on the basic nacelle model. A list of these configurations and their basic geometric 
parameters is given in table I. Three shrouds and six plugs were used and these are 
shown in the photographs in figure 4. Sketches giving dimensions and pressure-orifice 
locations on the shrouds and plugs are presented in figure 5. 
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The configuration designation has been coded as follows: 

(1) The first character, a number, indicates the boattail angle in degrees. 

(2) The second character, a letter, indicates the type of base used on the plug as 

follows: 

R no base 

F flat base 

S semitoroidal base 

(3) The third character, a number, indicates the amount of plug truncation in per- 
cent of full-length cut off. The 72.5-percent truncation is rounded to 73 percent. 

Figure 6 presents photographs of several of the configurations investigated. 

In strum entati on 

The gross thrust minus drag of the nozzle system was measured by a one- 
component, strain-gage thrust balance shown in figure 1. Static pressure was measured 
on the plug and afterbody surfaces at the locations shown in figure 5 by use of electronic 
pressure transducers. In addition, four total-pressure probes, as shown in figure 2, 
were used to measure the jet total pressure. Jet total temperature was measured at the 
location shown in figure 2 by a chromel-alumel thermocouple. Two turbine-type elec- 
tronic flowmeters located in the liquid hydrogen peroxide supply line were used to obtain 
liquid hydrogen peroxide mass flow. The electrical signals of all instruments were 
recorded by an automatic magnetic tape-recording system. 

Tests 

All configurations were investigated at static conditions (M = 0) and at free-stream 
Mach numbers from 0.50 to 1.30. However, some data were lost due to bad pressure 
orifices and excessive zero shifts of the thrust balance measurement. The angle of 
attack was held at a constant value of 0° during the entire investigation. The ratio of 
jet total pressure to free-stream static pressure was varied from 1.0 (jet off) to 
about 8.0. 


DATA ANALYSIS 


The data, recorded as electrical signals on magnetic tape, were used to compute 
standard force and pressure coefficients with an electronic computer. Pressure coef- 
ficients were referenced to free-stream dynamic pressure for Mach numbers greater 
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than 0 and to ambient pressure for static conditions (M = 0). All force coefficients were 
referenced to the maximum cross-sectional area of the model and free-stream dynamic 
pressure or ambient pressure. 

Approximately 15 frames of data were taken over a time period of about 0.75 sec- 
ond for each data point, and the average value was used for computations. Pressure 
forces on the boattail and cone plug were obtained by assigning to each pressure orifice 
an incremental area projected on a plane normal to the model axis and by numerically 
integrating the incremental forces. Model support interference is small for this support 
system and hence no correction for strut interference was made for this investigation. 
(See refs. 5 and 6.) 

The thrust balance measured the sum of the following axial forces: total momen- 
tum flux at the nozzle throat, plug pressure forces, external aerodynamic drag of the 
afterbody aft of station 104.39 cm, some internal tare pressure area forces in the nacelle, 
and friction forces. Thrust-minus-drag force for the nozzle was obtained by the use of 
the following equation: 


F - D - + (Pj - Poo)A-max + (^cyl^oAmax (1) 

The calculated skin -friction drag on the cylindrical portion of the shroud was not charged 
to the total drag of the nozzle system because the cylindrical portion of the shroud is not 
regarded as a part of the nozzle system. 

To obtain internal performance, which includes total momentum flux at the nozzle 
throat plus plug forces, the following equation was used. 


F - (F - D) + ^Cj-j^q^Amax + ^f^^oAmax 


( 2 ) 


For comparative purposes, the forces obtained by using equations (1) and (2) were 
nondimensionalized by dividing by the optimum isentropic thrust F^. This thrust is 
given by the equation 


F i 





( 3 ) 


and obtained by assuming low chamber fluid velocity and isentropic nozzle flow so that 
chamber temperature is equal to the stagnation temperature. Also, complete isentropic 
expansion of the exhaust flow is assumed so that p e equals p,^. Since jet total- 
pressure ratio and jet stagnation temperature are dependent on the efficiency of the 
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silver-screen catalyst bed used in the HgOg simulator, an average straight-line fairing 
of optimum isentropic thrust data from tunnel runs for all configurations was used to 
eliminate data scatter due to catalyst bed decay and error in throat areas. The maxi- 
mum deviation from this average line was approximately 4 percent; however, the majority 
of the data points were much closer to the average fairing. 

RESULTS AND DISCUSSION 
Afterbody Pressure Distributions 

Figure 7 presents the jet-off pressure distributions on the external surface of the 
three afterbodies for various Mach numbers. Pressure coefficients are shown for sev- 
eral configurations at the three radial rows of pressure orifices on each shroud (0°, 90°, 
and 180°). Several configurations are shown for the 8° and 20° boattail shrouds to indi- 
cate the effect of plug geometry on the jet-off pressure distributions. 

The good agreement of the pressure coefficients measured on the three radial rows 
of orifices on each afterbody indicates a symmetrical flow field around the model and also 
indicates that model support interference effects were small. Although no specific con- 
clusion can be made about the effect of plug geometry on afterbody pressure distributions 
from the data shown; the effect of plug geometry on the level of the pressure distributions 
is generally small or nonexistent, as might be expected. The pressure distributions 
shown in figure 7(a) for configuration 8-R-O were obtained with a shroud different from 
that used to obtain the pressure distributions for configurations 8-F-73 and 8-S-73. Due 
to warpage, the 8° boattail shroud was rebuilt early in the investigation and it is thought 
that small differences in rebuilding the second shroud produced the shift in the pressure- 
distribution levels shown in figure 7(a). Configuration 8-R-O was the only configuration 
investigated with the original 8° boattail shroud. 

The small dip in the pressure distribution on the cylindrical portion of the 20° 
boattail shroud as shown in figure 7(b) at supersonic Mach numbers was probably caused 
by some misalinement of the shroud at the point of attachment to the model 
( x /dmax = “0*?). 

All three shrouds exhibit the abrupt decrease in pressure over the initial boattail 
curvature characteristic of the external flow expansion around a corner or radius. Since 
disturbances can be propagated upstream through the boundary layer for shorter dis- 
tances at high speeds than at low speeds, the decrease in pressure on the boattail is more 
abrupt at supersonic speeds than at subsonic speeds. Pressure recovery is observed 
over most of the boattail length of the 8° and 20° boattail shrouds whereas external flow 
separation occurred on the 30° boattail shroud as shown in figure 7(c). After separation, 
the boattail pressures remained essentially constant. 
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Figure 8 shows the jet -on pressure distribution on the external surface of the three 
afterbodies for various Mach numbers and jet total-pressure ratios. Since the pressures 
measured around the afterbodies at 0°, 90°, and 180° locations agreed closely at the same 
longitudinal locations as shown in figure 7, the pressure coefficients are presented as 
averages in figure 8. The abrupt decrease in pressure over the external boattails and 
the external flow separation on the 30° boattail shroud shown for the jet-off pressure 
distributions of figure 7 are also exhibited for the jet-on pressure distributions. These 
pressure distributions are typical for boattail bodies. (See refs. 5 and 7.) 

Boattail pressures on all three shrouds decreased with an initial increase of jet 
total -pressure ratio. However, when jet total -pressure ratio is increased to higher 
values, the boattail pressures increased with increasing jet total pressure. The boat- 
tail pressure decreasing and then increasing with increasing jet total -pressure ratio was 
reported in reference 5 and is caused by the jet plume initially aspirating the boattail 
and base region at low jet total -pres sure ratios and then expanding into the free-stream 
flow at higher jet total -pres sure ratios, thus reducing the aspiration effects. Positive 
boattail pressure coefficients, and hence beneficial thrust terms, are obtained on the 8° 
boattail shroud at the higher values of jet total -pressure ratio whereas the pressure 
coefficients generally remained negative on the 20° and 30° boattail shrouds. 

Boattail Drag 

Boattail pressure drag, obtained by integrating the pressures over the axially pro- 
jected boattail area, is presented in coefficient form in figures 9 and 10. The variation 
of jet -off boattail drag coefficient with Mach number for the three shrouds investigated is 
shown in figure 9 for several configurations. The results of reference 8 indicate some 
model interference from the bow shock wave reflected from the tunnel wall in the region 
1.05 < M < 1.13. However, no interference from the bow shock wave is apparent for the 
data presented. A decrease of jet-off boattail pressure drag is observed for the 8° boat- 
tail shroud at M = 1.00. This decrease of drag is probably caused by some tunnel wall 
interference associated with the model nose-cone — cylinder juncture expansion field. 

This effect at M = 1.00 was observed in the Langley 16-foot transonic tunnel and 
reported in reference 8 for a model 152.4 cm in length and 15.24 cm in diameter. 

The three shrouds investigated have approximately the same projected boattail area 
and hence the jet-off drag is only a function of the expansion of the external flow around 
the boattail. For subsonic Mach numbers, the 8° boattail shroud exhibits the lowest 
pressure drag level, as might be expected, since the external flow is required to expand 
through a larger angle and to lower pressure levels for the 20° and 30° boattail shrouds. 
However, from M = 0.90 to M = 1.30, the lowest drag level is exhibited by the 30° 
boattail shroud. This low drag level is caused by separation of the external flow on the 
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30° boattail at these speeds, which prevents complete expansion around the boattail and 
avoids the attendant low pressures and high drag. The 8° boattail shroud exhibits 
approximately the same jet -off pressure drag level at Mach numbers above 1.1 as the 
30° boattail shroud. 

Considering the entire Mach number range of the investigation, the 8° boattail 
shroud generally exhibits the lowest jet-off boattail drag level whereas the 20° boattail 
shroud exhibits the highest boattail drag level. The amount of plug truncation and the 
plug base shape have little or no effect on the jet -off boattail drag levels. 

The effects of jet operation on boattail drag are shown in figure 10. Truncation of 
the cone plug up to 72.5 percent of full length had little or no effect on the jet -on boattail 
drag of the 8° boattail shroud. However, the 20° boattail shroud exhibits an increase in 
jet -on boattail drag when the cone plug was truncated from 72.5 to 100 percent of full 
length and a semitoroidal base was installed. Only one plug configuration was investi- 
gated with the 30° boattail and no conclusion can be made concerning the effect of plug 
geometry on the drag of this shroud. 

Jet operation initially increases boattail pressure drag on all three shrouds inves- 
tigated. All three shrouds also show some decrease of boattail drag in the high range of 
jet total-pressure ratios. The amount of drag decrease is dependent on boattail angle. 

Jet operation becomes favorable (i.e., the boattail drag at some jet total -pressure ratio 
is less than the jet-off value of boattail drag) on the 8° boattail shroud for jet total- 
pressure ratios larger than approximately 3.5 at all test Mach numbers. Jet operation 
is detrimental on the 20° boattail shroud throughout the range of jet total-pressure ratios 
for subsonic Mach numbers with the possible exception of the highest jet total-pressure 
ratio. At supersonic speeds, operation of the jet again has a favorable effect on the boat- 
tail drag of the 20° boattail shroud at the higher jet total -pres sure ratios. Jet operation 
had a detrimental effect on the boattail drag of the 30° boattail shroud for all Mach num- 
bers and jet total -pressure ratios of the investigation. 

A hysteresis effect is exhibited by the boattail drag of the 20° boattail shroud as jet 
total-pressure ratio is varied continuously (i.e., being increased and then decreased) at 
sonic and low supersonic speeds. This effect was reported in reference 5 and is attrib- 
uted to the axial movement of the external flow separation point on the boattail. Although 
separation also occurred on the 30° boattail shroud, the separation point apparently 
remained at a fixed point and hence no hysteresis effect is present on the 30° boattail 
shroud. 
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Cone -Plug Static -Pressure Distributions 

Figures 11, 12, and 13 present the cone-plug static-pressure distributions for sev- 
eral Mach numbers, jet total -pressure ratios, and configurations. These plug pressure 
distributions are typical of plug -type nozzles. (See refs. 5, 9, and 10.) The series of 
increases and decreases of local plug static pressure along the length of the cone plug, 
particularly at low Mach numbers, are a result of alternating regions of expansion and 
compression of the jet flow starting with a Prandtl-Meyer expansion fan originating at 
the boattail lip. The average plug base pressure, for truncated plugs, is shown in fig- 
ures 11, 12, and 13 as a solid symbol and generally has approximately the same level as 
the last pressure measured on the plug surface. 

The theoretical pressure distributions calculated by the method of characteristics 
as given in the appendix are compared with the experimental pressure distributions in 
figures 14, 15, and 16. The theoretical jet -plume shapes are also given in these figures 
for each case computed. The location of shock waves or regions of strong compression 
predicted by the method of characteristics is indicated in the figures by a dash-dot -dash 
line in the jet flow field. 

At regions of strong compression, deviation from the experimental distributions 
occurs because the flow becomes nonisentropic. Two such regions of strong compres- 
sion were predicted for configuration 8-R-0, as shown in figure 14(a). The location of 
the second strong compression may not be accurate because of errors introduced by the 
first region of strong compression. For the 72.5 -percent -truncated plugs and jet total- 
pressure ratios greater than 2.77 (fig. 16), the region of strong compression did not 
occur on the plug surface and hence the pressure distributions for these cases are not 
affected by nonisentropic flow. An abrupt deviation of the jet boundary at the intersec- 
tion of the boundary with a region of strong compression is also predicted by theory. 

Theory deviates from experimental data in the initial expansion region because no 
iterations were performed and the initial line was assumed to be straight. However, the 
main cause of disagreement between theory and experiment in this region is probably the 
assumption of flow direction perpendicular to the geometric throat which was inclined 15° 
from the normal (perpendicular to plug surface). From this assumption and two- 
dimensional calculations, the flow deflection angle at the initial line ^Mj = 1.02^ did not 
vary much from the value of -15° at the geometric throat. In actuality, for the 8° boat- 
tail shroud configurations, the flow deflection angle varied from -8° at the shroud lip to 
-15° at the plug surface along the geometric throat and also should have approximately 
the same variation along the Mj = 1.02 initial line. A comparison of the pressure dis- 
tributions for the 20° boattail shroud configuration (15° internal boattail angle) shown in 
figure 16(b) and the 8° boattail shroud configurations (fig. 16(a)) shows better agreement 
between theory and experiment for the 20° boattail shroud configuration. For the 20° 
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boattail shroud configuration, the actual flow deflection angle approaches the theoretical 
flow deflection distribution of -15° along the geometric throat. 

After the initial decrease in the theoretical static -pres sure distribution, a discon- 
tinuity of slope occurs and the static pressure on the plug tends to remain fairly con- 
stant. This region of the jet flow field consists of intersecting compression and expansion 
characteristic lines of opposite families whereas the preceding region of the jet flow field 
consisted almost entirely of expansion characteristic lines. 

The next region of the jet flow field consists almost entirely of compression char- 
acteristic lines and accounts for the abrupt increase in the static -pressure distribution 
shown by experiment and predicted by theory. In this region, strong compressions or 
shock waves begin forming and the theoretical pressures calculated downstream of this 
point have errors introduced by nonisentropic flow. 

Cone-Plug Base Pressure Distributions 

The radial distributions of cone -plug base pressure are presented in figures 17 
to 21 for several Mach numbers, jet total -pres sure ratios, and configurations. The pres- 
sures obtained on the base of the 100-percent -truncated plug configurations were sym- 
metrical about the plug center line. Therefore, the pressures at corresponding radial 
stations about the plug center line were averaged and are shown as such in figures 20 
and 21. Some asymmetric pressure distributions were obtained on the longer plugs and 
hence the base pressures measured on each side of the plug center line are shown for 
those configurations utilizing the 30-percent and 72. 5-percent -truncated plugs. 

Plug base pressure generally decreased at most conditions with increasing plug 
truncation. For example, a comparison of data for the 8° boattail flat -base configura- 
tions shown in figures 17, 18, and 20 shows the plug base region to be generally pres- 
surized for the 30 -percent -truncated plug, pressurized or aspirated depending on jet 
total -pres sure ratio for the 72.5 -per cent -truncated plug, and aspirated by the jet flow for 
the 100-percent -truncated plug. 

Variation of the shroud boattail angle and, hence, of the internal lip angle had little 
or no effect on the 72. 5-percent -truncated plug base pressures. However, a large effect 
of internal angle was observed for the 100-percent -truncated plug configurations at sub- 
sonic Mach numbers. This effect might be expected since large boattail and internal 
shroud angles direct both the external flow stream and the jet flow stream toward the 
plug base and thus tend to pressurize this region. Figure 21 indicates that the increase 
in base pressure of the 30° boattail over that of the 20° boattail is not as large as the 
increase in base pressure of the 20° boattail shroud over that of the 8° boattail. Boat- 
tail and internal shroud angle had little or no effect on plug base pressures at Mach 
numbers greater than 0.90. 
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The plug base shape had little effect on the plug base pressures, and no consistent 
trend is evident from the data shown. However, a comparison of the flat- and 
semitoroidal-base configurations indicates that the semitoroidal base established and 
maintained a vortex -ring type of flow more readily than the flat plug base. A vortex- 
ring type of flow is indicated by pressure peaks near the center and at the edges of the 
plug base. (See ref. 11.) Little or no vortex flow was established on the 30-percent- 
truncated plug. Vortex flow was maintained on the 72. 5-percent -truncated plug at most 
test conditions. Vortex flow was indicated on the 100-percent -truncated, flat -base plug 
only at the higher values of jet total -pres sure ratio and on the 100-percent -truncated, 
semitoroidal-base plug at all conditions at subsonic Mach numbers and at the higher 
values of jet total-pressure ratio at transonic and supersonic Mach numbers. 

The plug base pressures measured on plug configurations having some plug surface 
for jet flow expansion (truncation < 100 percent) increased or decreased depending on the 
value of jet total-pressure ratio. The plug base pressures are hence a function of the jet 
flow expansion along the length of the plug, as might be expected from the average plug 
base pressures shown in figures 12 and 13. The plug base pressures on the 100-percent- 
truncated plug configurations generally decreased with increasing jet total-pressure 
ratio. One exception is configuration 30-S-100 at subsonic speeds for which increasing- 
jet total -pres sure ratio increased the plug base pressures. 

Plug Thrust 

Figure 22 presents the variation of plug thrust coefficient with jet total -pres sure 
ratio and Mach number for all configurations investigated. Plug thrust was obtained by 
integrating the pressures over the axially projected plug area downstream of the annular 
throat. 


Nozzle Performance Characteristics 

Basic data .- Figure 23 presents the variation of the static thrust -minus-drag coef- 
ficient and the ideal static thrust coefficient with jet total -pressure ratio for all config- 
urations tested. 

Thrust -minus -drag coefficients, thrust coefficients, and ideal thrust coefficients 
for all configurations are shown in figure 24 as functions of jet total -pres sure ratio for 
the Mach numbers and jet total pressures of the investigation. It should be noted that the 
ideal -thrust-coefficient line is a faired, or averaged, line obtained from the ideal thrust 
data for all configurations investigated. The same faired line of the ideal thrust coef- 
ficient is used as the common reference for each configuration in figure 24. This pro- 
cedure merely averages any scatter in mass-flow and jet -temperature measurements. 
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Thrust coefficient was not obtained on several configurations due to loss of pressure data 
from which external drag is obtained. 

Thrust -minus -drag ratio .- The variation of thrust -minus -drag ratio with jet total - 
pressure ratio and amount of plug truncation is shown in figures 25 to 28 for several 
Mach numbers and configurations. The thrust -minus -drag ratio generally increases 
with increasing jet total-pressure ratio. Due to the increase of external drag with 
increasing Mach number, thrust -minus -drag ratio generally decreases with Mach number 
for the 0- and 30-percent -truncated plugs. Similar trends are shown for isentropic plug 
nozzles in reference 5. The variation of thrust-minus -drag performance with Mach num- 
ber for the 72.5- and 100-percent -truncated plugs was more erratic at transonic and low 
supersonic Mach numbers due to the larger effect of the plug base pressures. Thrust- 
minus-drag ratio decreases with increasing amount of plug truncation. Large losses 
were incurred by all configurations utilizing the 72.5- and 100-percent -truncated plugs. 
However, the thrust -minus -drag performance of the 30-percent -truncated plug configu- 
ration approaches the thrust -minus -drag performance of the full-length plug except at 
Mach numbers of 1.00 and 1.15. Generally, the 15° cone plug may be truncated as much 
as 30 percent of full length with little or no loss of thrust -minus -drag performance; fur- 
ther increase in plug truncation, however, incurs large losses on the nozzle thrust-minus- 
drag performance. 

Figures 29 to 32 present the effect of jet total -pres sure ratio and boattail angle on 
thrust-minus-drag ratio for several Mach numbers and configurations. At Mach numbers 
up to 1.15, and, in particular, for the 7 2. 5 -per cent -truncated plug configurations, the 
8° boattail configurations generally exhibit the best thrust -minus -drag performance char- 
acteristics due to lower external drag. The 100-percent -truncated, semitoroidal-base 
plug configurations shown in figure 32 exhibit trends which may be expected from exami- 
nation of the plug base pressures shown in figure 21. The 20° boattail gives the best 
performance due to higher plug base pressures at Mach numbers of 0 and 0.50. At Mach 
numbers above 0.50, the 8° boattail configuration generally gives the highest performance. 
The 30° boattail configuration is not shown at Mach numbers of 0.50 and 0.70 due to exces- 
sive zero shifts of the data. 

The effect of jet total-pressure ratio and plug base shape on thrust -minus -drag 
ratio is shown in figures 33 to 36 for several Mach numbers and configurations. A defi- 
nite trend of thrust -minus -drag ratio with plug base shape is not evident from the data 
shown. The semitoroidal base configurations usually exhibit the highest thrust-minus- 
drag performance for 72. 5-percent -truncated configurations and the lowest performance 
for the 100 -percent -truncated configurations. However, several exceptions are shown in 
figures 33 to 36. The effect of plug base shape appears to be dependent to a large extent 
on the flow field and the formation of the ring -type vortices; hence, determination of the 
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best plug base shape should be made in the flow field of the airplane model in which the 
nozzle system is to be installed. This requirement is not uncommon for some nozzle 
systems. (See ref. 12.) 

Thrust ratio .- The variation of thrust ratio or internal performance with jet total- 
pressure ratio is shown in figures 37 to 43 for several Mach numbers and configurations. 
The internal performance was obtained by adding the external pressure drag and skin 
friction to the thrust -minus -drag values. Values of thrust ratio are not shown for sev- 
eral configurations because pressure instrumentation was lost and hence no external 
pressure drag could be calculated. Thrust ratio is also not given at a Mach number of 0 
since, in this case, external pressure drag and skin friction equal zero and thrust ratio is 
equivalent to the thrust -minus-drag ratio shown previously. 

Internal performance generally follows the same trends as thrust -minus -drag per- 
formance since the thrust term constitutes the major terms in both parameters. Internal 
performance is generally at a higher level, however, since external drag is not charged 
to the internal performance. Internal performance generally increases with increasing 
jet total -pressure ratio and decreases with increasing Mach number. 

The effect of plug truncation on internal performance is shown in figures 37 and 38. 
Increasing the amount of plug truncation decreases the internal performance due to 
adverse pressures on the plug bases. The loss incurred by cutting the plug 30 percent 
of full length is small except at Mach numbers of 1.00 and 1.15. 

Figures 39 to 41 present the effect of boattail angle on internal performance. Inter- 
nal performance decreased with increasing boattail angle with the exception of the 
72. 5-percent -truncated flat-base plug configurations shown in figures 39(c) at M = 1.20 
and M = 1.30. 

The effect of plug base shape on internal performance is presented in figures 42 
and 43. Internal performance for the 100-percent -truncated configurations is not given 
because of the loss of pressure instrumentation for these configurations. The 
semitoroidal-base plug configurations exhibit the highest internal performance at Mach 
numbers up to 1.15 for the configurations shown. No conclusion can be made on the 
effect of plug base shape at Mach numbers of 1.20 and 1.30. 

Performance at a sche duled pressure ratio .- Figure 44 presents a jet total- 
pressure-ratio schedule for a typical turbojet engine for external Mach numbers from 
0.5 to 1.3. (See refs. 5 and 13.) The variation of nozzle internal performance and thrust- 
minus-drag ratio with Mach number for this typical turbojet total-pressure-ratio sched- 
ule is presented in figure 45. 

Figure 45(a) compares the performance of the full-length 15° cone plug with that of 
an isentropic plug reported in reference 5 and of a 10° cone plug nozzle reported in 
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reference 13. The isentropic plug nozzle was designed for the same conditions as the 
cone plug nozzle of the present investigation. The isentropic plug incurred much greater 
losses than the full-length cone plug. The low performance level of the isentropic plug 
is a result of high drag on the steep boattail angle required for isentropic plugs and of a 
large amount of overexpansion losses on the plug as evidenced by the low level of thrust 
ratio or internal performance in figure 45(a) (see ref. 1). The full-length 15° cone plug 
shows a higher level of performance as theorized in reference 10 because of the low 
value of boattail angle and less overexpansion penalty on the plug due to the external 
flow. The performance of the 15° cone plug nozzle is somewhat below the performance 
of the 10° (/3 = 5.5°) cone plug nozzle, except for Mach numbers slightly greater than 1.2, 
because of higher boattail drag and overexpansion losses on the plug. The nozzle per- 
formance of all configurations decreases with increasing subsonic Mach numbers and 
increases with increasing supersonic Mach numbers except for the 10° cone plug nozzle. 

The effect of plug truncation on nozzle performance for the typical turbojet total- 
pressure-ratio schedule is shown in figure 45(b) as a function of Mach number. Thrust 
ratio is not shown for configuration 8-F-100 because pressure instrumentation was lost 
for this configuration. Both thrust -minus-drag ratio and thrust ratio decrease with 
increasing amount of plug truncation. However, the loss incurred by truncating the plug 
30 percent of full length is generally less than 1 percent except in the transonic speed 
range. The 72.5- and 100 -percent -truncated plugs exhibit large losses as a result of 
plug truncation. 

Figure 45(c) presents the effect of boattail angle on nozzle performance for the 
typical turbojet total-pressure-ratio schedule. As a result of pressure instrumentation 
loss, thrust ratio is not shown for configuration 8-S-100, Cone plug nozzle performance 
generally decreases with increasing boattail angle. For Mach numbers less than 0.80, 
the 8° and 20° boattail shroud configurations exhibit approximately the same thrust- 
minus-drag performance. This is a result of the increase of plug base pressures shown 
in figure 21 for these Mach numbers due to increasing boattail angle. The resultant 
increase in plug thrust is greater than the increase of boattail drag due to a larger boat- 
tail angle, and the thrust -minus -drag performance of configuration 20-S-100 is hence 
increased to a level favorably comparable to that of configuration 8-S-100 (generally less 
than 1 percent loss). In the transonic speed range where the drag rise occurs, the 
increase in boattail drag becomes larger than the increase in plug thrust due to increased 
boattail angle, and thrust-minus-drag performance decreases with increasing boattail 
angle. Mixed results were obtained at a Mach number of 1.30. 
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CONCLUSIONS 


An investigation of the effect of plug truncation, boattail angle, and plug base shape 
on the performance of a 15° cone plug nozzle has been made. The nozzle was designed 
for a jet total-pressure ratio of 16.5 and a ratio of nozzle throat area to maximum cross- 
sectional area of 0.25. The investigation was conducted at Mach numbers from 0 to 1.30, 
jet total -pres sure ratios from 1.0 to approximately 8.0, and an angle of attack of 0°. The 
following trends are indicated: 

1. Nozzle performance loss as a result of a 30-percent truncation of the 15° cone 
plug was generally less than 1 percent whereas further truncation incurred large losses 
in nozzle performance. 

2. Cone -plug -nozzle performance generally decreased with increasing boattail 

angle. 

3. Although plug base shape had a large effect on nozzle performance, a definite 
trend of cone -plug-nozzle performance with plug base shape is not evident from this 
investigation. 

4. The 15° cone-plug-nozzle performance was higher than that of a comparable 
isentropic plug nozzle but somewhat lower than a 10° cone plug nozzle. 

5. Beneficial jet interference was observed at high values of jet total -pressure 
ratio on the 8° boattail shroud pressure; the boattail pressures remained negative at all 
conditions for the 20° and 30° boattail shrouds. 

6. Plug base pressures appear to be a strong function of the jet flow expansion 
process along the length of the plugs which were not fully truncated. For the 100- 
percent -truncated plugs, base pressure was a strong function of boattail angle at sub- 
sonic Mach numbers; base pressure and plug thrust increased with increasing boattail 
angle in the subsonic speed range. 

Langley Research Center, 

National Aeronautics and Space Administration, 

Langley Station, Hampton, Va., December 17, 1968, 

720-03-00-01-23. 
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METHOD - OF - CH ARAC T E RISTIC S COMPUTER PROGRAM 
FOR CONE PLUG NOZZLES 

The method of characteristics (refs. 14, 15, 16, 17, and 18) was used to calculate 
the theoretical plug pressure distributions presented in the body of this report. The 
computer program used for this purpose is reproduced as written in Fortran IV language. 
Also given are the program assumptions, restrictions, and inputs. Figure 46 is a sketch 
showing the initial and boundary conditions assumed for the cone plug nozzle. 

Assumptions 

The following conditions were assumed for the computer program used: 

1. Axisymmetric, supersonic, nonviscous, irrotational, isentropic, steady flow 

2. Perfect gas 

3. Jet flow is parallel everywhere along the geometric throat and inclined to the 

X-axis at an angle equal to the plug half -angle 

4. Two-dimensional flow at corner of shroud and along throat 

5. Throat or initial starting line is assumed downstream of geometric throat so 

that the initial Mach number is greater than 1.00 

6. No iterations 


Restrictions 

Use of the computer program is restricted by the following conditions: 

1. Cone plug nozzles only; no plug curvature 

2. No internal jet expansion 

3. No external flow 

Initial and Boundary Conditions 

A straight line inclined with the geometric throat at a Prandtl-Meyer turning angle 
corresponding to Mj = 1.02 was used for the initial conditions. This line is a charac- 
teristic line of the (0 - ju) family. Additional initial conditions were needed to define the 
shape of the characteristic line of the (0 + ijl) family. These initial conditions were 
obtained by generating a two-dimensional Prandtl-Meyer expansion fan from the boattail 
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lip. At the intersections of the characteristic line of the (8 + /x) family and the Prandtl- 
Meyer waves, the flow properties were calculated from two-dimensional isentropic rela- 
tions and, thus, the initial conditions along the (8 + ju) characteristic line were defined. 
The flow deflection angle 8 was allowed to vary lineally along the initial (8 - p) charac- 
teristic line from the value (approximately 15°) calculated at the first intersection of the 
(8 - p) and (8 + p) characteristic lines to the value on the plug surface (15°). 

Boundary conditions that must be met occur at the plug surface (flow deflection 
angle) and at the outer boundary at the exhaust stream (pressure). The ratio of specific 
heats used for all calculations was 1.266 (for products of decomposition of 90-percent 
hydrogen peroxide) and 1.0133 x 10^ N/m^ was used for atmospheric pressure. 

Characteristic Nets 

Figure 47 shows two examples of the characteristic net obtained by the method of 
characteristics for no external flow. The characteristic net and the initial and boundary 
conditions are presented in figure 47(a) for configuration 8-F-73 at a jet total -pressure 
ratio of 4.07. Figure 47(b) presents the characteristic net calculated for configuration 
8-F-30 at a jet total-pressure ratio of 4.03. The transition from isentropic flow to non- 
isentropic flow is clearly predicted by the theory for this case. A shock wave is formed 
when compression waves of the same family converge and intersect. (See ref. 15.) The 
shaded region shown in figure 47(b) is a region of strong compression of the jet flow as 
predicted by theory and which may form a shock wave in the actual flow case. These 
characteristic lines converged and intersected for several cases with lower jet total- 
pressure ratios calculated for this investigation and might for this case if a smaller 
mesh size were selected. An abrupt change of jet boundary at the strong compression or 
shock is also predicted. 


Input Cards 


The 
form with 


inputs required for the program are as follows, 
six decimal places unless otherwise specified. 


All input is in floating-point 


Card 

Columns 

Contents 

1 

1 to 5 

II = an integer 

If II = 0, Mach number along initial (6 + ju) 
characteristic line (see fig. 46) read in 
as data 



If II = 1, Mach number calculated by 
program 


20 



APPENDIX 


Card 

Columns 

1 

6 

to 

14 


15 

to 

23 


24 

to 

32 


33 

to 

41 


42 

to 

50 


51 

to 

59 


60 

to 

64 


65 

to 

69 

2 

1 

to 

9 

3 

1 

to 

9 



10 to 18 


19 to 27 


28 to 36 

4 

1 to 10 


11 to 20 


21 to 30 


31 to 40 


41 to 50 


4(a) optional (Card 1 must have II = 0) 
1 to 9 

10 to 18 


Contents 

GAMMA = ratio of specific heats 
PR = jet total-pressure ratio p^. jjPoo 
YO = radius of plug base, inches 

FL = length of plug from shroud lip to plug 
base, inches 

YS = radius of shroud at throat, inches 

PGANRD = plug half -angle, radians 

ICONFG = an integer signifying configuration 
number 

ICASE = an integer signifying case number 

THRTAN = throat angle (negative), radians 

XMACH(1,1) = Mach number along initial 

( 9 - ji) characteristic line 
(see fig. 46) 

XMACH(2,1) = Mach number as above 

XMACH(3,1) = Mach number as above 

XMACH(4,1) = Mach number as above, located 
on plug surface 

TL = length of reference plug from shroud lip 
to plug tip, inches 

PLUGAN = plug half -angle, degrees 

TRUN = percent of plug truncation 

BETA = shroud internal angle, degrees 

XS = axial coordinate of shroud lip (generally 
taken as zero), inches 

XMACH(1, 1) = Mach number along initial 

(0 + i u.) characteristic line 
(see fig. 46) 

XMACH(1,2) = Mach number as above 
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Card Columns Contents 

4(a) optional (Card 1 must have n = 0) 


19 to 27 

XMACH(1,3) = 

28 to 36 

XMACH(1,4) = 

37 to 45 

XMACH(1,5) = 

46 to 54 

XMACH(1,6) = 


Mach number as above 

Mach number as above 

Mach number as above 

Mach number as above, located 
on jet boundary 


Computer Program 

The method -of-characteristics computer program was written and debugged by use 
of the English gravitational system of units. The program as written in Fortran IV lan- 
guage is as follows: 



oooooo ooooooooo 


METHOD OF CHARACTERISTICS - MAIN DECK 
CONE PLUG NOZZLE PRESSURE DISTRIBUTION 

SUPERSONIC, STEADY, ISENTROPIC, INVISCID, AXI -SYMMETR! C FLOW 

INITIAL LINE AND IDENTIFICATION SECTION 


DIMENSION X(48,45),Y(48,45),THETA(48,45), XMU( 48, 45 ),W(48,45),XL(48 
1,45),XM(48,45) , XRAT 10(48,45) ,XMACH( 48, 45 ) , PRAT I0(48,45),C(7,7),D(7 
2,7) 


SUBPROGRAM TO DETERMINE INITIAL VALUES FOR MAIN DECK 
SUBDECK TO METHOD OF CHARACTERISTICS 


REAQ(5, 1) I I , GAMMA, PR, YO,F L, YS, PGANRD, ICONFG, I CASE,THRTAN 
1 FORMAT (I5,6F9.6, 2I5/F9.6) 

READ (5, 920) (XMACH (1,1), 1=1, 4) 

920 FORMAT (4F 9 . 6) 

READ(5,4)TL, P LU GAN, TRUN, BETA, XS 
4 FORMAT ( 5F10 . 6) 

DPRATO = PR 

IONFIG = ICONFG 

IF (I I.EQ.DGO TO 902 

READ ( 5,903) (XMACHU, J), J = l,6) 

903 FORMAT ( 6F9 .6 ) 

I F ( I I . EQ. 0)G0 TO 904 

902 XMACH (1,6) =SQRT ( 2 . / ( GAMMA- 1 . ) * ( PR* * ( ( GAMMA- 1 . ) /GAMMA ) - 1 . ) ) 
XMACHd, 1) = 1.00 
DELTA=( XMACH ( l,6)-1.0)/5.0 
DO 905 J=2, 6 


to 

co 
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to 


K=J-1 

905 XMACH ( 1, J ) =XMACH (1,K)+DELTA 
XMACH(1,1)=1.02 

904 A=-.5*SQRT(( GAMMA+ 1 . ) / ( GAMMA- 1 . ) ) 

B= ( 1 . - GAMMA* * 2 ) / 2 . 

00- 906 = 1. ..K 

C( 1/ J ) = XMACH ( 1 / J)**2/ ( 1 . + ( GAMMA- 1 . )/ 2 . * XMACH ( 1, J ) ** 2 ) 

D( 1, J ) =1 ./ C( 1, J ) 

THETA(1,J)=A*(ARS!N(B*C(1,U) +GAMMA) -1 . 570788)-. 5*(ARSIN( GAMMA- 2 . *D 
1(1, J))+1.570788)+THRTAN 

XMU ( 1, J ) =ARS I N ( 1 . / XMACH ( 1, J ) ) 

906 W(l, J)=SQRT((GAMMA+1. )/2. *XMACH( 1, J )** 2/ ( 1 . + ( GAMMA-1 . )/2. *XMACH(1, 
1J ) **2 ) ) 

XX=(Y0+FL*TAN( PGANRD ) -YS) / (TAN( PGANRD )+TAN(THETA(l, 1 )-XMU( 1,1))) 
YY=Y0+(FL-XX)*TAN( PGANRD ) 

ALPHA=- (THETA (1,1 ) -XMU (1,1)+1. 570788) 

H=(YS-YY)/ COS (ALPHA) 

0ELH=H/4. 

DO 907 1=1,4 
X I = I 

DELY=XI *DELH*C0S( ALPHA) 

Y ( 1,1) =YS-DELY 

907 X ( I , 1 ) =-DELY*TAN (ALPHA) 

DO 8 J=2, 6 

K= J- 1 

X(1,J)=(Y(1,K)-YS-X(1,K)*TAN(THETA(1, K)+XMU ( 1, K) ) )/ (TAN (THETA ( 1,J ) 
1-XMU (1,J))-TAN(THETA(1,K ) +XMU ( 1, K ) ) ) 

8 Y(l, J)-Y(1,K) + (X(1, J)-X(1,K))*TAN(THETA(1,K)+XMU(1,K)) 

DE LTHE= ( -PGANRD -THETA ( 1, 1 ) ) / 3 . 

DO 909 1=2,4 
J = ! -1 

THETA( !,1)=THETA(J,1)+DELTHE 

W( !,1)=SQRT((GAMMA+1. )/2.*XMACH( I , 1 )**2/ ( 1. + (GAMMA-1 . )/ 2 . *XMACH( I , 
11 ) * *2 ) ) 

909 XMU ( 1,1) =ARS I N ( 1 . / XMACH ( 1,1)) 
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PRINT 910, GAMMA, PR/ F L, PGANRD, THRTAN, I CONFG, I CASE 

910 FORMAT ( 1H1, 2X, 5HGAMMA, 5X, 10HPRES RAT I 0, 4X, 6HPLUG L,3X,12HPLUG ANG 
1RAD,12H THO ANG RAD, 3X, 6HC0NF I G, 7X, 4HCASE/ /5F12 . 7, I 7, I 12 ) 

PRINT 911, (X(1,J),Y(1,J),W(1,J),THETA(1, J),XMU(1,J),J=1,6) 

911 FORMAT (1H0,14H INITIAL I LINE///2X, 1HI ,4X, 1HJ, 7X,1HX, 11X,1HY, 11X,1H 
1W,9X,3HTHETA,9X,2HMU//2X,1H1,4X,1H1,2X,5F12.7/2X,1H1,4X1H2,2X,5F12 
2.7/2X,lHl,4X,lH3,2X,5F12.7/2X,lHl,4X,lH4,2X,5F12.7/2X,lHl,4X,lH5,2 
3X,5F12.7/2X,1H1,4X,1H6,2X,5F12.7) 

PRINT 912, (X(I,1),Y(I,1),W(I,1), THETA (1,1), XMU( 1,1), 1=1, 4) 

912 FORMAT (1H0, 14HINITIAL J LI NE///2X, 1H I ,4X, 1HJ, 7X, 1HX, 11X, 1HY, 11X, 1H 
1W,9X,5HTHETA,9X,2HMU//2X,1H1,4X,1H1,2X,5F12.7/2X,1H2,4X1H1,2X,5F12 
2.7/2X,lH3,4X,lHl,2X,5F12.7/2X,lH4,4X, 1H1, 2X, 5F12. 7) 

PRINT 5, I ON FIG, DP RATO, TRUN, BETA, GAMMA, PLUGAN, FL,TL 
5 FORMAT ( 1H1, 7X, 64HCALCULAT I ON OF LOCAL PLUG PRESSURES BY METHOD OF 
1 CHARACTER I STICS////2X,6HCONFIG,2X,10HPRES RAT 1 0, 3X, 4HTRUN, 6X, 4HBET 
2A, 6X, 5HGAMMA, 2X, 10HPLUG ANGLE, 2X, 7HFRACT L,3X,5HREF L// 3X, I 4, 4X, F6 
3. 2, 5X, F 6 . 2, 4 X , F 6 . 2,4X,F6. 3, 3X, F6 . 2, 4 X, F 7 . 3, 3X, F 7 . 3 ) 

PRINT 10, ( X C 1 ,1),Y( I, 1), THETA ( I , 1 ) , XMU ( I , 1 ) , W ( I , 1 ) , I =1, 4 ) 

10 FORMAT (1H0,22X,35HINITIAL J LINE CONDITIONS I =1, 2, 3, 4//5X, 1HX, 10X, 
11HY, bX, 9HTHETA RAD,2X,9HM ANG RAD, 5 X, 4HV/ A*// ( 5F 11 . 7 ) ) 


CALCULATION OF PRANQTL-MEYER EXPANSION NAVES AND PLUG INTERSECTION 
PRANDTL-MEYER WAVE SECTION 


DO 100 J=2, 6 
L = J-l 

N = J + 3 
M = N-l 
DO 200 1=2, M 
K = 1-1 

X( I , J) = (Y ( I , L)-Y(K, J)+X(K, J ) *TAN (THETA (K, J)-XMU(K, J))-X( I, L)*TAN(T 
1HETA( I, L)+XMU( I , L ) ) ) / ( TAN ( THETA ( K, J ) -XMU ( K, U ) ) -TAN (THETA ( I , L)+XMU ( 
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2 I , L) ) ) 

Y(!,J)=Y(I,L)+(X(I,J)-X(I, L))*TAN(THETA( I , L)+XMU( I , L) ) 

XL( I , L) =S I N( XMU ( I , L) ) *S I N (THETA ( I , L ) ) *TAN ( XMU ( I , L ) ) / COS (THETA(I/L) 

1 + XMU ( ! , L ) ) 

XM( K, U ) =S I N ( XiMU ( K, J ) ) *S I N ( THETA ( K, J ) ) *TAN ( XMU ( K, J ) ) / COS ( THETA ( K, J ) 
1-XMU (K, J ) ) 

XTHETA = (1.-W( I ,L)/W(K,J)-(THETA( I ,L)-THETA(K,J) ) *TAN (XMU( K, J ) ) + (X( 
1I,J)-X(K,J) )*XM(K, J )/Y(K, J)-XL( l,L)/Y(l,L)*(X(!,J)-X(l,L))*W(l,L)/ 
2W ( K, J ) ) / ( TAN ( XMU ( K , J ) ) +W ( I , L ) / W ( K, J ) *TAN ( XMU ( I , L ) ) ) 

THETA( I , J )=THETA( ! , D+XTHETA 

DW=W( I , L) * ( TAN (XMU ( I , L) )*XTHETA+XL( l / L)/Y(! / L)*(X(l,J)-X(f / L))) 
W(! / J)-W(I / L)+DW 

XMACH ( ! ,J)=W( ! / J)*(1./(SQRT((GAMMA+1. )/2.-( (GAMMA-1. )/ 2. )*W( l,J)** 
12))) 

XII = II 

I F ( ( XMACH (I, J). LT. 1.00). AND. (XI I. GT. 2. 5)) GO TO 611 
I F( XMACH ( I , J ) . LT. 1.00) GO TO 506 
ARG=1. 0/(SQRT( XMACH ( I , J)**2-l. 0) ) 

I F ( XMACH ( I , J ) . GE . 1.00) GO TO 200 

506 PRINT 507 

507 FORMAT (1H1////37H MACH NUMBER IS SUBSONIC -- NOT VALID///49H READ 
1MACH ( l y 1 ) THROUGH MACH(1,6) INTO SUBPROGRAM// 74H MACH ( 1, 1 ) = 1.02, 

2 MACH (1,2)= 1.09 AND MACH(1,6) DEPENDS ON PRESSURE RATIO//65H MACH 
3 1 NE IS RECALCULATING USING THESE VALUES AND ONES SHOWN BELOW) 

XMACH ( 1 , 6 ) =S QRT ( 2 . / ( GAMMA- 1 . ) * ( PR * * ( ( GAMMA- 1 . ) /GAMMA ) - 1 . ) ) 

XMACH ( 1, 1 ) =1 . 02 
XMACH ( 1, 2 ) =1 . 09 

BBZ= (XMACH (1,6 ) -XMACH ( 1,2))/ 4.0 
XMACH (1,3) = XMACH ( 1, 2 ) +BBZ 
XMACH (1,4) -XMACH (1,3)+BBZ 
XMACH (1,5) =XMACH ( 1, 4 ) +BBZ 
PRINT 601, ( XMACH ( 1, K),K=1,6) 

601 FORMAT ( 1H /// 60H MACH(1,1) MACH(1,2) MACH(1,3) MACH(1,4) MACH(1,5) 
1 MACH ( 1, 6 ) / / 6F 10 . 6 ) 

I 1=3 
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GO TO 904 

611 PRINT 612 

612 FORMAT ( 1H ///81H MACHINE HAS RECOMPUTED USING MACH VALUES ABOVE AN 
IQ MACH-NUM&E& IS ST LU SUBSONlCy/lOH NOJ VAL \AUJ35\\ MACH I NE I S JL 
2RED AND IS THEREFORE CONCLUDING OPERATION) 

GO TO 11 

200 XMU(l,d)=ATAN(ARG) 

PRINT 2, (X( I ,d ), Y( I ,d ),THETA( I ,d ) ,XMU( l,d),W(l,d),l=l,M) 

2 FORMAT (1H0,12X, 54HCALCU LATED FIELD POINTS FOR ONE J LINE AND 1=1,2 
1, 3. . .M//5X, 1HX, 10X, 1HY, 6X, 9HTHETA RAD,2X,9HM ANG RAD, 5X, 4HV/A*// ( 5 
2F11.7)) 

PRINT 3, J,M 

3 FORMAT ( 1H0, 4X, 1HJ, 8X, 1HM/ / I 6, 4X, 15) 

PLUG INTERSECTION OF PRANDTL-MEYER WAVES AND STATIC PRESSURES 

X(N,d)=(YO+X(M,d )* TAN (THETA (M, J )-XMU(M, J ) ) +F L*TAN ( PGANRD )-Y (M, J ) ) / 
1 ( T AN ( PGAN RD ) + TAN ( TH E T A ( M, d ) - XMU ( M , d ) ) ) 

Y(N,d)=YO+(FL-X(N,d))*TAN (PGANRD) 

THETA (N, J )=- PGANRD 

XTHETA=THETA(N, J ) -THETA (M, J) 

DW=W (M, d ) * ( XM(M, d )/ Y (M,d )* (X (N, d ) -X (M, d ) ) -TAN ( XMU (M, J ) ) *XTHETA ) 

W(N, d ) =W (M, d ) +DW 

XMACH (N,d )=W(N,d)*(l./( SQRT ( ( GAMMA+ 1 . ) / 2 . - ( (GAMMA-1 . ) / 2 . ) *W( N, d ) ** 
12))) 

ARG=1.0/(SQRT(XMACH(N,d)**2-1.0)) 

XMU(N,d)=ATAN(ARG) 

XRAT IO(N,d)=X(N,d)/TL 

PRAT I 0(N,d) = ( 1.0- (GAMMA-1. )/ (GAMMA+1 . )*W(N,d )**2)** (GAMMA/ (GAMMA-1 
1.)) 

PRINT 7,N,d,X(N,d),Y(N,d),THETA(N,d), XMU (N,d),W(N,d), XMACH (N,d),XR 
lATIO(N,d), PRAT I 0 ( N, d ) 

7 FORMAT ( 1H0, 23X, 32HCALCULATED POINT ON PLUG SURFACE// /2X,lHI,4X,lHd 
1 / 7X,1HX,9X,1HY,4X, 9HTHETA RAD,2X,9HM ANG RAD, 3X, 4HV/A*// I 3, I 5 , F 12 . 
26,4F10.6///1X,10HLOCAL MACH, 2X, 7HX/ REF L,3X,7HP/P TOT//3F10.6) 


DO 

-3 
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to 

00 


100 CONTINUE 
C 
C 

C INITIAL CALCULATION ON PRESSURE BOUNDARY -- 2 - D I M . SECTION 

C 

X ( 1, 7)=(Y(1,6)-YS+XS*TAN(THETA(1,6))-X(1,6)*TAN(THETA(1,6 ) +XMU (1,6 
1) ) )/ (TAN (THETAC 1/ 6) ) -TAN (THETA ( 1,6 ) +XMU ( 1,6 ) ) ) 
Y(1,7)=X(1,7)*TAN(THETA(1,6))-XS*TAN(THETA(1,6))+YS 
WW=W(1,6) 

XL(1,6)=SI N (XMU (l,6))*SIN(THETA(l,6))*TAN(XMU(l,6))/COS(THETA(l,6) 
1+XMU(1 # 6)) 

XTHETA= ( WW-W ( 1, 6 ) ) / ( W( 1, 6 ) *TAN ( XMU ( 1, 6 ) ) ) -X L( 1, 6 ) * ( X ( 1, 7 ) -X( 1, 6 ) )/ 

1(Y(1/G)*TAN (XMU (1,6))) 
v'l(l,7)=W(l / 6) 

THETA ( 1/ 7 ) =THETA( 1, 6 )+XTHETA 
XMU (1,7) =XMU (1,6) 

PRINT 14, X ( 1, 7 ) , Y ( 1, 7 ) ,THETA( 1, 7 ) , XMU( 1, 7 ) ,W( 1, 7 ) 

14 FORMAT (1H1,23X,34HINITIAL POINT ON PRESSURE BOUNDARY/// 2X, 1H I , 4X, 1 
1HJ,7X,1HX,9X,1HY,4X,9HTHETA RAD,2X,9HM AMG RAD, 3X, 4HV/A*/ / 2X, 1H1, 4 
2X,1H7,F12,6,4F10,6) 


PRESSURE BOUNDARY CALCULATION SECTION 


DO 300 J=7, 45 
L= J-l 
N=J + 3 
M-N-l 
I l =0-5 

l F ( J . LE . 7 ) GO TO 12 

I I = J-6 
J J = J- 7 

X(II,J) = (Y(I I, L)-Y(JJ,L) + X(JJ,L)*TAN(THETA(JJ, L ) ) — X ( I I , L)*TAN (THET 
1A( I I , L ) + XMU ( I I, L)))/(TAN(THETA(JJ,L))-TAN(THETA(I 1 , L) +XMU ( I 1 , L) ) ) 


APPENDIX 



o o o o o 


Y ( ! I , J ) =X ( I l,J)*TAN(THETA(JJ,L))-X(JJ,L)*TAN(THETA(JvJ,L))+Y(JJ,L) 

X L( I I , L)=SIN(XMU( I I, L) )*S IN (THETA ( I I , L) )*TAN(XMU( I I / L))/COS(THETA( 
II l,L)+XMU(l I , L) ) 

W(l I ,<J )=W( JJ, L) 

XTHETA=(W( I I, J)-W( I I, L))/(W( I I , L)*TAN (XMU( I I , L) ) ) -XL( I I , L)* (X ( I I , J 
1 )-X C I l / L))/(Y(l I , L)*TAN(XMU( I l,L))) 

THETA( I I , J )=THETA( I I, D + XTHETA 
XMU(I I, J)=XMU(JJ,L) 

PRINT 13, I I , J , X ( I I , J ) , Y ( I I, J) /THETA ( I I , J ) , XMU ( . I I , J ) ,W< I I , J ) 

13 FORMAT (1H0,21X,37HCALCU LATED POINT ON PRESSURE BOUNDARY/ //2X, 1HI ,4 
IX, 1HJ/ 7X,1HX,9X, 1HY, 4X,9HTHETA RAD,2X,9HM ANG RAD, 3X, 4HV/A*// 13,15 
2/F12.6/4F10. 6) 


CALCULATION OF FIELD POINTS SECTION 


I I = J- 5 

12 DO 400 1 = 1 I/M 

K= I - 1 

X( l,J) = (Y( I, L)-Y(K,J) + X(K, J ) *TAN ( THETA (K/J) -XMU (K/ J ) )-X( I / L)*TAN (T 
1HETA( I , L)+XMU( I , L) ) )/ (TAN( THETA ( K,J )-XMU( K, d )) -TAN (THETA ( I , L)+XMU( 
21, L))) 

Y(I,J)=Y(I,L)+(X( I , J ) -X ( !,L))*TAN(THETA(!,L)+XMU(I,L)) 

XL( I , L) =S I N(XMU ( l,L))*SIN(THETA( I , L) ) * TAN (XMU (l,L))/COS( THETA ( I , L ) 
1+XMU ( I , L)) 

XM(K, J ) =S I N (XMU ( K, J ) ) *S I N (THETA ( K,J ) ) *TAN (XMU( K, J ) )/ COS (THETA ( K,J ) 
1-XMU ( K, J ) ) 

XTHETA=(1.-W( I , L)/W(K, J)- (THETA ( I , L)-THETA(K, J ) )*TAN(XMU(K, J ) )+(X( 
II ,d)-X(K,J))*XM(K,J)/Y(K,J)-XL( l,L)/Y< l,L)*(X( l,J)-X( l,L))*W( l,L)/ 
2W ( K, J ) ) / ( TAN ( XMU ( K, J ) ) +W ( I , L ) /W ( K,'J ) *TAN ( XMU ( I , L) ) ) 

THETA( I , J )=THETA( I , D+XTHETA 

DW=W( I , L)*(TAN(XMU( I , L) )*XTHETA+XL( I , L)/Y( I ,L)*(X( l,J)-X( I, L))) 

H(I,J)-W(I,L)+DW 

XMACH ( l,J)=W( I , J )* ( 1 . / (SORT ( ( GAMMA+ 1 . ) / 2 . - ( ( GAMMA- 1 . ) /2 . )*W( I , J )** 


to 

CD 
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CO 

o 


12)3) . 

ARG=1.0/(SQRT(XMACH( I ,d )**2-1.0) ) 

400 XMU(I ,d)=ATAN(ARG) 

PRINT 2, (X( I ,d ),Y( l,d),THETA( I , J ) , XMU ( l , J ) ,W ( I , J ) , I = I I ,M ) 
PRINT 9,0,11 ,M 

9 FORMAT (1H0,4X,1HJ,8X,2HI I , 8X, 1HM// I 6, 4X, I 5, 4X, 15) 


FIELD CHARACTERISTICS INTERSECTION ON PLUG SECTION 


X(N,d)=(YO+X(M,d)*TAN(THETA(M, d ) -XMU (M,d )) +F L*TAN ( PGANRD ) -Y(M, d )) / 
1 ( TAN ( PGANRD ) +TAN ( THETA ( M, J ) -XMU (M, J ) ) ) 

IF(X(N,d).GT.FL)GO TO 11 
Y(N,J)-YO+(FL-X(N,J))*TAN( PGANRD) 

THETA (N,d)=- PGANRD 
XTHETA=THETA(N, J) -THETA <M,d) 

DW*W(M,d)*(XM(M,d)/Y(M,d)*(X(N,d)-X(M, J ) ) -TAN ( XMU (M, d))*XTHETA) 

W (N, J ) =W (M, d ) +DW 

XMACH (N,d)=W(N,d)*(l./(SQRT((GAMMA+l. )/ 2. -((GAMMA-1. )/2. )*W(N,d)** 
12 ))) 

ARG=1. 0/ (SQRT ( XMACH (N, d ) **2-1 . 0) ) 

XMU(N, d ) =ATAN ( ARG ) 

XRAT IO(N,d)=X(N,d)/TL 

PRAT IO(N,d)=(1.0-( GAMMA- 1 . ) / ( GAMMA+ l.)*W(N,d)**2)**( GAMMA/ (GAMMA-1 
1 . )) 

PRINT 7,N,d,X(N,d),Y(N,d),THETA(N,d), XMU (N,d),W(N,d), XMACH (N, d ) , XR 
1AT 1 0 ( .N, d ) , PRAT I 0(N, d ) 

300 CONTINUE 
11 STOP 
END 
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Figure 1.- Sketch of hydrogen peroxide jet simulator. All dimensions are in centimeters. 
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Figure 2.- Sketch of cone plug nozzle installed on pylon nacelle. All dimensions are in centimeters. 
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Figure 3.- Photograph of cone-plug-nozzle pylon nacelle installed in test section. L-66-4083 
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Figure 4.- Photographs of cone plugs and shrouds. 
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Figure 5.- Sketch of cone plugs and shrouds giving dimensions and orifice locations. All dimensions are in centimeters. 
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Figure 6.- Photographs of several cone-plug-nozzle configurations. 
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(a) p = 8°. 

Figure 8.- Effect of jet total-pressure ratio on average boattail pressure distributions for various Mach numbers. 
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(b) (5 = 20°. Concluded. 
Figure 8.- Continued. 
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Figure 8.- Concluded. 
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Boattail pressure drag coefficient. 
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Local plug pressure ratio, Ppiug/Pf j 



«) M = 1.00. 
Figure 11.- Continued. 
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Local plug pressure ratio, P p | ug /Pt 
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(a) P tj /P„ = 2.87. 

Figure 14.- Comparison of theoretical and experimental plug static-pressure distributions on configuration 8-R-O. Dash-dot-dash lines repre- 
sent location of shock waves or regions of strong compression; symbols with tails indicate row 180° orifices. 
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Figure 15.- Comparison of theoretical and experimental plug static-pressure distributions on configuration 8-F-30. Dash-dot-dash lines repre- 
sent location of shock waves or regions of strong compression; symbols with flags indicate row 180 orifices. 
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(a) Configuration 8-F-73. 

Figure 16.- Comparison of theoretical and experimental plug static-pressure distributions on configurations 8-F-73 and 20-F-73. Dash-dot-dash lines represent 
location of shock waves or regions of strong compression; symbols with flags indicate row 180° orifices. 
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(b) Configuration 20-F-73. 
Figure 16.- Concluded. 
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Configuration 8-R-O 



(a) Configuration 8-R-O. 

Figure 22.- Variation of plug thrust coefficient with jet total-pressure ratio and Mach number. Tailed symbols indicate decreasing 

pressure ratios. (Note shifting scales.) 
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Configuration 8-S-IOO 



(h) Configuration 8-S-100. 
Figure 22.- Continued. 
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Static thrust-minus-drag coefficient, C^ F _ D ^, and ideal static thrust coefficient, 



(a) p = 8°; flat base. 

Figure 23.- Variation of static thrust-minus-drag coefficient and ideal static thrust coefficient with jet total -pressure ratio 
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(f) Configuration 8-S-100. 
Figure 24.- Continued. 
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(i) Configuration 20-F-100. 
Figure 24.- Continued. 
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(k) Configuration 30-S-100. 
Figure 24- Concluded. 
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Jet total- pressure ratio, p. 


(b) M = 0.80 and M = 0.90. 
Figure 27.- Continued. 
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(c)M = LOO and M = 1.15. 


Figure 27.- Continued. 
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(d) M = 1.20 and M = 1.30. 


Figure 27.- Concluded. 
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Configuration 20-S — 100 
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Jet total -pressure ratio, pj. j/Pqo 
(b) M = 0.80 and M = 0.90. 


Figure 28.- Continued. 
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Figure 30.- Effect of jet total- pressure ratio and boattail angle on thrust-minus-drag ratio. 
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(b) M = 0.80 and M = 0.90. 
Figure 30.- Continued. 
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Jet total-pressure ratio, p t j/Poo 


(a) M = 0 to M = 0.70. 

Figure 32.- Effect of jet total- pressure ratio and boattail angle on thrust-minus-drag ratio. 100-percent truncation; semitoroidal plug base. 
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(d) M = 1.20 and M = 1.30. 
Figure 33.- Concluded. 
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Figure 34.- Continued. 





Configuration 20 F 73 
— Configuration 20— S-73 
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Jet total -pressure ratio, p f j/Poo 

(c) M = LOO and M = 1.15. 


Figure 34.- Continued. 
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Configuration 20-F — IOO 




(a) M = 0 to M = 0.70. 

Figure 36.- Effect of jet total-pressure ratio and plug base shape on thrust-minus-drag ratio. 20° boattail angle; 100-percent truncation. 
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Configuration 20-F — 100 


Configuration 20~S — 100 
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Jet total-pressure ratio, p t j/Poo 


(c) M = LOO and M = 1.15. 


Figure 36.- Continued. 
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Configuration 8-F-73 

Conf iguration 8-S-73 
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Jet total -pressure ratio, pj. j/Poq 


(c) M = 1.20 and M = 1.30. 
Figure 42.- Concluded. 
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Jet total-pressure ratio, j/Pqq 


(a) M = 0.50 to M = 0.80. 

Figure 43.- Effect of jet total-pressure ratio and plug base shape on thrust ratio. 20° boattai! angle; 72.5-percent truncation. 
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Conf iguration 20-F-73 

Configuration 20— S-73 



(c) M = 1.20 and M = 1.30. 
Figure 43.- Concluded. 
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Thrust ratio, F/F\ Thrust-minus-drag ratio, (F-D)/Fj 


Configuration 8-S-I00 
Configuration 20-S-100 
Configuration 30-S-I00 
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(c) Effect of boattail angle. 
Figure 45.- Concluded. 





Figure 46.- Sketch showing boundary and initial conditions required for method-of-characteristics computer program. 
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Configuration 8-F-73 



(a) Configuration 8-F-73. 

Figure 47.- Characteristic net for cone plug nozzles. 
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Configuration 8-F-30 
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(b) Configuration 8-F-30. 
Figure 47.- Concluded. 
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